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Ganadería, protección del hábitat y diversidad de abejas y avispas 
en el desierto del Monte Central 


H RESUMEN. El principal objetivo de las reservas es prevenir o mitigar los 
impactos humanos sobre los ecosistemas naturales. Es importante evaluar cuán bien 
las reservas alcanzan este objetivo. Evaluamos si la protección del hábitat que brinda 
la Reserva de la Biósfera de Nacufián (Monte Central, Argentina) resulta en cambios 
detectables en la estructura del hábitat, y en la riqueza y la composición de especies 
de abejas y avispas. Realizamos muestreos con trampas bandeja y observaciones 
de visitantes florales en seis pares de sitios dentro y fuera de la reserva. Nuestros 
resultados sugieren que los treinta y cinco años de exclusión del ganado vacuno en 
Nacufián han tenido efectos detectables sobre la estructura del hábitat. Sin embargo, 
estos cambios en el hábitat se tradujeron sólo en efectos parciales y conflictivos sobre 
la riqueza de himenópteros, y no tuvieron efectos detectables sobre la composición de 
himenópteros. Nuestro estudio debería repetirse en el futuro, con un mayor esfuerzo 
de muestreo y a lo largo de varios años antes que estos resultados puedan ser aplicados 
como guía de decisiones de manejo. 


PALABRAS CLAVE. Abejas. Avispas. Ganadería. Composición y riqueza de especies. 
Protección del hábitat. 


a ABSTRACT. Reserves aim mainly at preventing or mitigating human impacts 
on natural ecosystems. It is important to assess how well reserves meet this goal. We 
evaluated whether habitat protection offered by Ñacuñán Biosphere Reserve (Central 
Monte desert, Argentina) results in detectable changes in habitat structure and the 
species richness and composition of bees and wasps. We conducted pan trap sampling 
and flower visitor observations in six pairs of protected and unprotected sites. Our 
results suggest that thirty five years of cattle exclusion in Nacufián have had detectable 
effects on habitat structure. However, these changes in habitat structure translated 
only into partial and conflicting effects on hymenopteran richness, and did not have 
detectable effects on hymenopteran composition. Our study should be repeated with a 
greater sampling effort and throughout multiple years before our results can be applied 
to guide management decisions. 


KEY WORDS. Bees. Wasps. Cattle raising. Species composition and richness. Habitat 
protection. 


INTRODUCTION of human activities such as agriculture, 

logging and livestock grazing in natural 

Nature reserves are usually created to ecosystems. Thus, it is important to evaluate 
prevent or mitigate the widespread impact how well reserves meet their goal. However, 
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although reserves may help the conservation 
of target habitats or species of special 
concern, they are probably too few, small 
and isolated to represent a sufficient solution 
for the preservation of biodiversity, making it 
necessary to incorporate human-dominated 
landscapes into conservation strategies. 
Under proper management practices, 
natural habitats under livestock grazing and 
other types of human impacts may indeed 
be adequate for an integrative strategy to 
biodiversity conservation (Gómez-Aparicio 
et al., 2004; Steffan-Dewenter et al., 2006; 
Kremen et al., 2007; Winfree et al., 2007). 

Human-dominated ecosystems may 
be important for conservation in their own 
right. In the specific case of cattle raising, 
there is much controversy around its impact 
on natural ecosystems, particularly in arid 
regions (Fleischner, 1994; Noss, 1994; 
Wuerthner, 1994; Brown & McDonald, 
1995; Bestelmeyer & Wiens, 2001; Tabeni 
& Ojeda, 2005). Most livestock species 
are exotic to the ecosystems in which they 
graze, and thus can be seen as a widespread 
human-caused, chronic disturbance. 
Alternatively, because many of the original 
native grazers-under whose influence native 
species have evolved-have become much 
rarer than in the past or have gone extinct, 
the presence of livestock may resemble 
more the original disturbance regime than 
their absence. Regardless of our stand in 
this debate, it is again important to know 
how human activities and their exclusion 
influence the structure and dynamics of 
ecological systems. 

Bees and wasps (Hymenoptera) comprise 
an important group of pollinators and 
potential enemies of insect pests and are 
particularly abundant and diverse in arid 
regions (Michener, 1979, 2000). They 
perform important community processes 
and ecosystem services as pollen vectors, 
predators and parasites. Variation in the 
diversity of bees and wasps can be related 
to changes in the structure and abundance 
of floral and nesting resources (Gess & 
Gess, 1991; Samejima et al., 2004; Lassau 
& Hochuli, 2005; Potts et a/., 2005; Loyola 
& Martins, 2008). Browsing and trampling 


by livestock can contribute to changes 
in vegetation and habitat structure, and a 
reduction of grazing pressure through habitat 
protection can affect the diversity of bees and 
wasps (Kruess € Tscharntke, 2002; Sjódin, 
2007). 

Nacufián Biosphere Reserve (Mendoza, 
Argentina) provides an excellent opportunity 
to evaluate the effectiveness of habitat 
protection from livestock grazing for 
biodiversity conservation. Located atthe heart 
of the Central Monte Desert, Nacufián is a 
12,000 ha rectangular reserve with 35 years 
of exclusion of human activities, surrounded 
by active cattle ranches (Boshoven et 
al., 2001). Here we ask whether habitat 
protection offered by the reserve results in 
detectable changes in habitat structure, and 
whether such changes, if observed, affect in 
turn the species richness and composition of 
bees and wasps. 


MATERIAL AND METHODS 
Study area and sites 


Our study region lies within the Central 
Monte desert of Argentina (Cabrera, 1971; 
Roig-Juñent et al., 2001). The climate is 
semi-arid and most rainfall occurs in the 
spring and summer (October-March). We 
worked in Nacufián Biosphere Reserve 
(34°02’ S, 67°58’ W) and surrounding areas. 
Average annual rainfall in the region is 280 
mm (Estrella et al., 2001). Grazing has been 
excluded from the reserve since 1972, but has 
continued outside (Boshoven et al., 2001). 
Although we have no good estimates of cattle 
burden in these ranches, we have repeatedly 
seen cattle and their signs of activity (e.g., 
feces, tracks) in our non-reserve study sites 
(see below). Logging of native species has 
been banned in the whole Province for 
many decades. The algarrobo woodland 
is the most representative community, in 
which Prosopis flexuosa is accompanied 
by the small tree Geoffrea decorticans and 
several shrub (Larrea divaricata, Condalia 
microphylla, Capparis atamisquea, Lycium 
spp.) and grass species (Pappophorum spp., 
Trichloris crinita, Digitaria californica; Roig, 
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1971). Trees and most shrubs in the area 
are insect-pollinated (D. P. Vázquez and V. 
Aschero, personal observation). 

We worked in six pairs of protected and 
unprotected sites, each pair consisting of two 
1-ha plots, one inside and one outside the 
reserve, at least 1 km away from the fence 
marking its limit. Distance between paired 
plots was smaller than distance within pairs. 


Habitat structure: plant cover and 
composition, floral abundance and soil 
compactness 


To characterize the effects of cattle on 
habitat structure we estimated overall plant 
cover and composition, floral abundance 
and soil compactness. The latter is relevant to 
determine habitat structure because many of 
the hymenopteran species in our study nest 
in the ground, and soil compactness could 
thus affect their nesting success. Vegetation 
cover is relevant because it may affect soil 
humidity and temperature (again influencing 
ground-nesting insects) and it may affect 
nesting sites for wood-nesting insects (e.g., 
greater cover could mean greater availability 
of nesting sites). Plant composition could 
also influence nesting site availability for 
some species (wood-nesting species). Finally, 
flower abundance is a measure of resource 
availability for flower-visiting insects (all 
bees and some wasps). Although these four 
components of habitat structure are arguably 
phenomenological, in the sense that we 
are not measuring the nesting and feeding 
requirements of each hymenopteran species 
in our study, this is the best we can do in 
a community-wide study involving many 
species. 

In each 1-ha plot we ran two transects 
crossing each other at the center of the plot. In 
each transect we estimated percent cover for 
each plant species and bare ground by direct 
observation in 20 m x 20 m quadrats every 
five meters (i.e., a total of 39 quadrats per 
plot, because the two transects overlapped 
in the center). We also counted the number 
of flowering individuals per quadrat of 
each animal-pollinated plant species as an 
estimate of overall floral abundance. Floral 


abundance was estimated periodicaly, at the 
same dates in which flower visitor sampling 
was conducted (see below). 

To assess soil compactness, we measured 
soil penetration resistance using a penetrometer 
device. Penetration resistance was defined 
as the number of hits needed for the device 
to penetrate 5 cm belowground. Penetration 
resistance was measured at each plot at twenty 
sampling points regularly spaced along 4 
equidistant parallel transects. 


Pan trap sampling 


We conducted pan trap samplings of 
hymenopterans in four different occassions 
during the spring of 2005 (late October, early 
and late November, and mid December). 
On each sampling date, nine pan traps 
were placed in each 1-ha plot to collect 
bees and wasps. Traps were placed at 
regular distances along the two transects 
used for plant sampling, and were exposed 
approximately from 8 am until 6 pm, when 
trapped specimens were collected from 
the traps. Each trap unit consisted of three 
small blue, white and yellow plastic trays, 
filled with water with a drop of detergent as 
surfactant. 


Flower visitor sampling 


We sampled insects visiting flowers of 
the most common entomophilous plant 
species every 15-20 days throughout the 
2005 flowering season (October-December). 
Paired sites were sampled simultaneously by 
different observers in the morning (10:00- 
13:00 h) and the afternoon (15:00-18:00 h). 
Observers were switched between sites in 
each sampling period to avoid observer bias. 
At each plot, the observer used an aspirator, a 
sweep net or a killing jar with ethyl acetate to 
sampleallinsectsvisiting during 5 mina group 
of flowers or inflorescences of the study plant 
species: Capparis atamisquea (Capparaceae), 
Condalia microphylla (Rhamnaceae), Larrea 
divaricata and L. cuneifolia (Zygophilaceae), 
and Prosopis flexuosa (Fabaceae); these five 
species represented the bulk of the insect- 
pollinated flowers available in Nacuñán 
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during the study period. Individual plants 
were selected haphazardly for sampling, 
attempting to sample as many flowering 
individuals as possible in different locations 
within each plot, so as to minimize spatial 
autocorrelation among samples. 


Data analyses 


To evaluate if floral abundance and soil 
compactness differed between reserve and 
non-reserve sites we used generalized linear 
mixed model analysis, with “protection” 
as a fixed factor with two levels (reserve 
and non-reserve), “site” as a random factor 
with six levels (sites 1-6), a Poisson error 
distribution and a log link function. Analyses 
were conducted using the Imer function in 
the Ime4 package of R statistical software, 
version 2.4.0 (http://www. r-project.org). We 
used Akaike's information criterion (AIC) 
to compare the fit of mixed models with 
and without the interaction between fixed 
and random factors. The model with the 
interaction term was considered superior 
only if its associated AIC value was at least 
two units lower than the AIC associated to 
the model without interaction (Burnham & 
Anderson, 2002). We also checked model 
deviance for overdispersion, using the 
“estimated scale” parameter, ®, returned 
by Imer (i.e., the model residual deviance 
over the residual degrees of freedom), 
which indicates departure from the assumed 
Poisson distribution of errors (Crawley, 
2005). A model is overdispersed if d>1, and 
underdispersed if «1. When for a given 
model we found a substantial departure 
from one in the value of ®, we re-ran the 
analysis using the quasi-Poisson distribution 
(by setting family=quasipoisson in the Imer 
function call). 

Cover was compared between reserve 
and non-reserve sites with a paired t-test; this 
design is appropriate because there was only 
one cover estimate per plot. 

We conducted sample-based rarefaction 
to evaluate if species richness differed 
between reserve and non-reserve sites, re- 
scaling richness estimates to individuals 
(Gotelli 8 Colwell, 2001). To compare 


species composition we used non-metric 
multidimensional scaling (MDS), and 
evaluated the significance of site ordination 
with Mantel permutation tests (Legendre 
& Legendre, 1998; Vazquez & Simberloff, 
2003), reporting the Mantel statistic for the 
original data (rg) and its 95% permutation 
confidence interval (rp). Rarefaction analyses 
and Mantel permutations were conducted 
with ad-hoc functions written in R (http:// 
www.r-project.org) and Octave (http://www. 
octave.org), respectively, available from the 
authors upon request. MDS analyses were 
conducted using the isoMDS in the MASS 
package of R. 


RESULTS 


Soil hardness tended to be greater outside 
than inside the reserve (Fig. 1). The AIC for 
the model without interaction was 99.52, 
whereas that for the model with interaction 
was 99.8. Thus, the inclusion ofthe interaction 
term does not improve model fit, and we can 
conclude that the effect of the interaction is 
not statistically significant. The main effect 
of management in the generalized mixed 
model without interaction was statistically 
significant (z=7.003, P<0.0001). The 
estimated scale for this model was ®=0.59, 
thus indicating underdispersion. To account 
for this departure from the assumption that 
@®=1, we ran the model again with the 
quasi-Poisson distribution; this change in 
the distribution did not affect parameter 
estimates and their associated significance, 
suggesting that the effect of management on 
soil hardness was indeed significant. 

Floral abundance did not differ significantly 
between reserve and non-reserve sites (mixed 
model with interaction: z=0.731, P=0.465, 
@=1.01). Plant cover was significantly greater 
inside than outside the reserve (paired t-test: 
t=3.26, P<0.05; Fig. 2). Site ordination in 
terms of plant composition separated reserve 
and non-reserve sites, although the trend 
was non-significant, as the Mantel statistic 
calculated for the original data falls within 
the permutation confidence interval (Fig. 3). 

Eighty-six species of Hymenoptera were 
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Fig. 1. Soil hardness in each study site. Site pairs 1-6 are shown from left to right. Grey boxes, reserve 
sites; white boxes, grazed sites outside the reserve. 
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Fig. 3. Reduced-space plot of first two coordinates resulting from nonmetric multidimensional scaling (MDS) 
of plant species composition. Black symbols, reserve sites; white symbols, non-reserve sites. Symbol shape 
indicates site pair. Permutation test to evaluate statistical significance of reserve protection on composition: 
Mantel statistic for observed data: r9=-0.13; permutation confidence interval: r52[-0.24,0.24]. 
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Fig. 4. Comparison of species richness between reserve (black circles) and non-reserve (white circles) 
sites. Error bars indicate confidence intervals of rarefied richness for the most abundant sample in each 
pair (see methods). (a) Hymenoptera from pan traps. (b) Flower visitors. 
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Fig. 5. Reduced-space plot of first two coordinates resulting from nonmetric multidimensional scaling 
(MDS) of insect species composition. (a) Hymenoptera from pan traps. (b) Flower visiting Hymenoptera. 
Black symbols, reserve sites; white symbols, non-reserve sites. Symbol shape indicates site pair. 
Permutation test to evaluate statistical significance of reserve protection on composition: pan traps: 
rg70.07; rp=[-0.24,0.24]; flower visitors: 0.03; rj=[-0.25,0.24]. 
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collected in our study sites (Appendices | and 
II). Rarefied richness of Hymenoptera caught 
in pan traps inside the reserve was equal or 
greater than outside the reserve (Fig. 4a). In 
contrast, rarefied richness of flower visitors in 
reserve sites was equal or lower than outside 
the reserve, except for site 1, where richness 
was greater inside the reserve (Fig. 4b). 
Species composition of Hymenoptera did 
not show a consistent, statistically significant 
pattern between reserve and non-reserve 
sites, for both Hymenoptera caught in pan 
traps and flower visitors (Fig. 5). 


DISCUSSION 


Our results suggest that thirty five years 
of cattle exclusion in Nacufián Biosphere 
Reserve have had detectable effects on habitat 
structure. Soil inside the reserve appears to 
be less compact and with greater plant cover 
than in neighboring non-reservesites. Because 
many species of bees and wasps nest in the 
ground, these differences in soil compactness 
and cover can indeed be interpreted as 
ecologically meaningful habitat differences. 
There was also a trend towards differing plant 
composition which, from the perspective of 
bees and wasps, could also mean differences 
in habitat quality and structure, but not in 
flower abundance. 

Changes in habitat structure translated 
into partial and conflicting effects on 
hymenopteran richness and did not have 
detectable effects on  hymenopteran 
composition. Whereas pan trap results suggest 
that the reserve has led to a modest increase 
of hymenopteran richness inside the reserve, 
the flower visitor results suggest the opposite 
trend. It is reasonable to expect differences 
in the species composition of samples from 
multiple sampling methods, given that each 
has its own bias and works best for specific 
groups of Hymenoptera (Kearns & Inouye, 
1993; Southwood & Henderson, 2000). The 
latter is likely to be true for our two sampling 
methods; in fact, only 4696 of Hymenoptera 
visiting flowers were also caught in pan traps. 
However, there is no reason to expect different 
trends in the responses of these different sets 


of species to a particular environmental 
effect, unless the sampling bias correlates 
with some sort of ecological or life history 
traits that can influence species' responses. 
For example, if species trapped with one 
technique used a different set of resources 
than those trapped with the other technique, 
and if resources were differently affected by 
habitat protection, then different responses of 
the groups associated to these two methods 
would be expected. The limited data available 
on the natural history of the Hymenoptera of 
the Central Monte do not allow us to assess 
the validity of these speculations. 

Our results should be interpreted with 
cautionifusedtoguidemanagementdecisions. 
First, we have data for only one year, which 
may not be representative of what happens 
in the long term, given the great seasonal and 
interannual variability of ecological systems. 
Second, even though we sampled repeatedly 
throughout the spring, our sampling effort 
may not have been sufficiently large to detect 
many rare species, making our richness 
and composition estimates incomplete. In 
fact, lack of completeness of the richness 
estimates is suggested by the rather steep, 
non-saturating rarefaction curves (not shown). 
Despite these limitations, it is at least possible 
that our results are indeed a reflection of 
reality and that the reserve creation and its 
associated habitat changes have had little 
effect on the hymenopteran fauna. Only 
future studies, conducted over multiple years 
and with greater sampling effort, will allow 
this question to be answered. 
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Appendix I. Hymenopteran taxa recorded in pan traps. 














Family Morphospecies name Unprotected Site Protected Site 

1 2 3 4 5 6 1 2 3 4 5 6 
Bees 
Andrenidae Acamptopoeum sp. 1 1 
Andrenidae Andrenidae sp. 1 1 
Andrenidae Calliopsis sp. 1 1 
Andrenidae Callonychium (Callonychium) sp. 1 30 42 167 92 23 2 52 30 48 47 15 13 
Andrenidae Callonychium (Micronychium) sp. 1 1 3 4 1 7 2 4 
Andrenidae Psaenythia phillanthoides 1 2 
Andrenidae Psaenythia tricolor 23 43 2 5 5 18 4 2 210 2 2 
Andrenidae Psaenythia sp. 1 1 1 1 2 1 1 
Andrenidae Spinoliella sp. 1 13 16 3 2 9 1 1 6 6 10 3 17 
Apidae Alloscirtetica arrhenica 431 112 508 89 143 77 139 217 418 90 169 98 
Apidae Alloscirtetica brethesi 4 5 5 1 2 5 13 10 1 
Apidae Apis mellifera 1 1 1 
Apidae Diadasia patagonica 2 
Apidae Doeringiella sp. 1 1 1 
Apidae Isepeolus octopunctatus 1 3 2 2 6 3 4 5 
Apidae Isepeolus vachali 
Apidae Isepeolus sp. 1 1 
Apidae Leiopodus trochantericus 1 3 1 
Apidae Leptometriella separata 1 7 4 1 2 1 3 2 6 
Apidae Melectoides cockereli 1 
Apidae Rhinepeolus rufiventris 1 1 1 2 
Apidae Svastra bombylans 3 2 3 4 1 
Apidae Svastrides zebra 1 6 1 1 
Apidae Xylocopa ordinaria 1 
Colletidae Colletes sp. 1 1 
Colletidae Leioproctus (Perditomorpha) stilborhinus 1 2 
Colletidae Leioproctus sp. 1 5 1 3 1 1 1 
Colletidae Leioproctus sp. 2 1 1 
Halictidae Augochloropsis sp. 1 1 1 1 2 4 
Halictidae Caenohalictus sp. 1 1 2 1 d 1 
Halictidae Dialictus sp. 1 1 2 31 69 25 19 16 15 1 
Halictidae Dialictus sp. 2 175 361 472 880 581 1075 480 404 70 200 735 146 
Halictidae Pseudoagapostemon singularis 3 4 2 1 
Megachilidae Anthidium sp. 1 3 2 
Megachilidae Anthidium sp. 2 1 
Megachilidae Megachile sp. 1 1 
Megachilidae Megachile sp. 2 1 2 2 
Megachilidae Pseudocentron sp. 1 2 1 1 2 
Wasps 
Braconidae Braconidae sp. 1 1 
Chrysididae Chrysididae sp. 1 1 1 1 
Crabronidae Crabronidae sp. 1 1 1 
Crabronidae Crabronidae sp. 2 1 
Crabronidae Crabronidae sp. 3 1 
Crabronidae Solierella sp. 1 1 1 
Crabronidae Solierella sp. 2 1 3 1 2 5 1 3 4 2 6 
Ichneumonidae Anomaloninae sp. 1 3 1 2 2 4 3 6 
Pompilidae Ageniella sp. 1 1 
Pompilidae Ageniella sp. 2 1 
Pompilidae Pompilinae sp. 1 3 1 2 1 1 
Pompilidae Pompilinae sp. 2 3 1 
Sphecidae Ammophila sp. 1 1 1 1 2 
Sphecidae Prionyx pumilio 1 1 
Sphecidae Tachysphex sp. 1 1 1 1 
Sphecidae Trypoxylon sp. 1 1 2 
Tiphiidae Calchaquila albinervis 1 1 1 2 1 2 2 
Tiphiidae Tiphia sp. 1 9 
Vespidae Hypalastoroides sp. 1 1 1 2 1 2 8 1 3 1 1 1 
Vespidae Vespidae sp. 1 1 1 
Vespidae Zeta canaliculata 2 1 


continúa 
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Appendix II. Hymenopteran taxa recorded visiting flowers. 








Family Morphospecies name Unprotected Site Protected Site 
1 2 3 4 5 6 1 2 3 4 

Bees 

Andrenidae Augochloropsis sp. 1 1 1 1 

Andrenidae Psaenythia phillanthoides 1 2 1 2 

Andrenidae Psaenythia tricolor 1 

Apidae Alloscirtetica arrhenica 1 1 

Apidae Alloscirtetica brethesi 1 

Apidae Apis mellifera 1 5 1 

Apidae Brachynomadinae sp. 1 1 

Apidae Centris sp. 1 1 

Apidae Centris brethesi 

Apidae Doeringiella sp. 1 1 

Apidae Rhinepeolus rufiventris 1 

Colletidae Caupolicana ruficollis 1 

Colletidae Colletes sp. 1 2 3 2 4 1 1 6 2 4 

Colletidae Colletes bicolor 1 

Colletidae Leioproctus (Kylopasiphae) pruinosus 1 1 1 7 

Colletidae Leioproctus (Perditomorpha) stilborhinus 1 2 5 1 35 

Colletidae Leioproctus (Perditomorpha) sp. 1 1 

Colletidae Londopria sp. 1 

Halictidae Dialictus sp. 2 2 4 1 1 1 1 

Megachilidae Anthidium sp. 1 4^ 03 2 2 

Megachilidae Epanthidium sp. 1 1 

Megachilidae Megachile sp. 1 1 1 1 2 

Megachilidae Megachile sp. 2 1 

Megachilidae Megachile sp. 3 1 

Wasps 

Crabronidae Podagritus sp. 1 

Crabronidae Rubrica nasuta 1 

Eumenidae Eumeninae sp. 1 1 1 

Eumenidae Hypalastoroides sp. 1 1 

Gasteruptiidae Gasteruption sp. 1 1 

Mutillidae Mutillidae sp. 1 1 6 1 2 1 3 1 

Pompillidae Anopleus sp. 1 1 

Pompillidae Pompillidae sp. 1 1 

Scoliidae Campsomeris sp. 1 1 

Sphecoidea Larridae sp. 1 1 

Thiphiidae Thynninae sp. 1 3 5 1 1 1 

Thiphiidae Tiphia sp. 2 1 

Vespidae Monobia sp. 1 1 


Vespidae Zethus sp. 1 1 





